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A proposal for a scalable, solid-state implementation of a 
quantum computer is presented. Qubits are fluorine nuclear 
spins in a solid crystal of fluorapatite [Ca5F(P04)3] with reso- 
nant frequencies separated by a large field gradient. Quantum 
logic is accomplished using nuclear-nuclear dipolar couplings 
with decoupling and selective recoupling RF pulse sequences. 
Magnetic resonance force microscopy is used for readout. As 
many as 300 qubits can be implemented in the realistic lab- 
oratory extremes of T = 10 mK and Bo = 20 T with the 
existing sensitivity of force microscopy. 

PACS numbers: 03.67.Lx, 76.60.Pc, 07.79.Pk 



An increasing number of theoretical developments have 
recently motivated the construction of a quantum com- 
puter with a large number of qubits. The potential for 
efficient simulation of quantum systems J |], the discov- 
eries of algorithms for fast factorization p | and database 
searching p| , and the development of quantum error cor- 
rection B have helped to spur a large number of pro- 
posals for experimental implementations of a quantum 
computer. However, any physical implementation of a 
quantum computer must battle the fact that well-isolated 
quantum systems are difficult to couple and to mea- 
sure, whereas the introduction of necessary couplings and 
probes leads to the devastating effects of decoherence. 

To date, the most successful experimental realization 
of a multi-qubit quantum computer (or at least a sim- 
ulation thereof) has been in room-temperature solution 
NMR |5|. Here, the spin-states of molecular nuclei in a 
solution are well isolated, as demonstrated by long ther- 
mal relaxation times (Ti) of many seconds. The nuclei 
of each molecule are weakly coupled by scalar couplings. 
Measurement without substantial decoherence is made 
possible by the large ensemble of ~10^^ uncoupled, iden- 
tical molecules in the solution. There is much debate, 
however, as to whether solution NMR currently does or 
ever will exhibit the signatures of truly quantum infor- 
mation processing g]. The presence of a large, mixed en- 
semble renders quantum entanglement and wave- function 
collapse of individual spins unmeasurable. Certain ma- 
nipulations allow the processing of an "effective pure 
state" in which the ensemble behaves nearly identically 
to a single quantum system pM]. However, it has been 
shown that the evolution of effective pure states of exist- 
ing solution NMR computers can be described without 
any entanglement, and either more qubits or much larger 



nuclear polarizations will be needed before entanglement 
is demonstrable Q . Even if we neglect the debate of what 
is "quantum" and what is not, the principal limitation 
due to the small polarizations in solution NMR quantum 
computation is that the usable signal decreases exponen- 
tially with the number of qubits, leaving this method 
unlikely to exceed the 10 qubit level without extensive 
modification H. 

In this letter, we propose an implementation using 
solid-state NMR. In a solid crystal there exists the poten- 
tial for polarization of the nuclear spins by simple cooling 
of the lattice to extremely low temperatures. The ther- 
mal relaxation time Ti can be dramatically lengthened 
when moving from a solution to an insulating, nonmag- 
netic crystal, and the nuclear-nuclear coupling can be 
faster. Hence, a solid crystal has great potential for quan- 
tum computation in terms of timescales Effl and scala- 
bility. 
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FIG. 1. The fluorapatite crystal, of length 10 fim and 
width roughly 1 pm, is mounted near the end of a cantilever. 
A dysprosium micromagnet with dimensions W = 10 ^m, 
H = 4 /xm, D = 400 /im is aligned parallel to the crystal 
planes. A field gradient of 1.4 T/^m is produced at the crys- 
tal, a distance of 2.07 ^ra above the magnet. The insert shows 
the F sites of the crystal structure. The different shades in- 
dicate different resonant frequencies. 

These advantages, however, come at a cost: several 
difficult and related problems arise in implementing a 
solid-crystal NMR computer. This letter, although not 
the first proposal for quantum computation in a crystal 
lattice [Q, provides new approaches to these problems. 
One problem is the complicated network of dipolar cou- 
plings in a solid, which must be suppressed to slow T2 



decoherence but selectively retained for logic. We ad- 
dress this problem with both "selective averaging" radio 
frequency (RF) pulse sequences and a well-chosen crys- 
tal structure. Another problem arises from the need to 
distinguish and detect nuclei in the periodic lattice of a 
crystal; for this we establish a very large, static, one- 
dimensional magnetic field gradient with a microfabri- 
cated, high-magnetization ferromagnet [|l2| . Only a small 
ensemble of nuclei may fit into the area over which the 
gradient field is homogeneous, and therefore a more sen- 
sitive means than standard inductive pick-up is needed 
to measure the nuclear spin states. We propose using 
magnetic resonance force microscopy (MRFM). 

A schematic for the solid-crystal quantum computer 
we propose appears in Fig. |l|. The quantum computer is 
an ensemble of N one-dimensional chains of n spin-1/2 
nuclei. Due to the field gradient, the resonant frequencies 
of the nuclear spins within a chain are different. Hence 
the secular component of the dipolar Hamiltonian which 
couples the zth spin to the jth spin within the Tnth chain 
is written f[l3| 
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where 7 is the gyromagnetic ratio {2n x 40 MHz/T for 
^®F), a is the distance between spins in a chain, and (j) 
is the angle between the chain of spins and the large 
applied magnetic field, which is taken to be in the z- 
direction. Other terms of the dipolar Hamiltonian av- 
erage to zero on a timescale of l/Awij or faster, where 
Awij — 7a|Vi3f I is the separation of the ith and jth 
resonant frequencies caused by the field gradient ^ Bf^ 
between them. The Hamiltonian of Eq. (^ may be 
"switched ofP' by applying a periodic succession of nar- 
row band tt pulses at, for instance, the ith resonant fre- 
quency p3[ . Simultaneous decoupling of more than two 
qubits may be accomplished by timing the selective tt 
pulses according to the entries of an appropriately sized 
Hadamard matrix; a pair of qubits may be selectively 
recoupled in order to implement two-bit gates |lj] . 

Non-disturbing measurement is possible in this scheme 
because, orthogonal to the chain direction, a nucleus of 
resonant frequency tOi has a large plane of copies with 
equal frequency. The coupling between copies has a dif- 
ferent form; the nth and ?Tith nucleus of identical reso- 
nant frequency LOi are coupled by the Hamiltonian ||l3|| 
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where A„j„ is the distance between the two nuclei in units 
of a, and Omn is the angle between the vector which con- 
nects them and the direction of the applied field. This 
coupling between copies must be "switched off' to pre- 
vent interference between ensemble members. Decou- 
pling these spins may be done by any number of pulse 



sequences which have been developed over the past 40 
years in solid-state NMR, the simplest of which is the 
WAHUHA pulse sequence. Such sequences can reduce 
dipolar broadening by more than three orders of magni- 
tude |l5). 



Fortunately, the different forms of Ti 
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allow constant decoupling of copies without adverse ef- 
fect on the manipulations of couplings within individual 
computers. The qubit-qubit coupling of Eq. (|l|) may 
in general be rescaled and rotated by the broadband 
WAHUHA-type sequence, but it may still be controlled 
for logic gates with suitably phased -k pulses. The broad- 
band WAHUHA-type pulse sequence which decouples all 
copies is unaffected by the addition of selective 7r-pulses, 
insofar as those 7r-pulses are sufficiently short in com- 
parison to the pulse spacing. This criterion is difficult 
to meet, since the tt pulses must last a time of about 
f/Awi.i+i to be frequency selective. A simple example 
to illustrate the decoupling scheme is shown in Fig. H. 
This example is limited by the adverse effect of the finite 
TT pulse-length on the WAHUHA decoupling, but longer 
and more sophisticated sequences can be designed with 
standard techniques to compensate for this and other lim- 
iting effects | |l5[ . 
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FIG. 2. An example RF pulse sequence. The first line 
(WHH) shows broadband decoupling of qubit copies; the nar- 
row, shaded pulses are appropriately phased 7r/2 pulses [ |l5[ , 
made sufliciently narrow (broadband) to affect all nuclei in 
the crystal. The subsequent three lines show selective pulses 
to decouple qubits; these wider pulses are sufliciently long 
(narrowband) to affect only individual planes. Single spin ro- 
tations, denoted by </>, may be performed in between cycles. 
Multi-bit gates such as controUed-NOT may be performed by 
combining single spin rotations with selective recoupling. For 
example, qubits i and i + 1 may be recoupled by adding a n 
pulse at the time noted by the dotted line |L4| . 

The only problematic dipolar couplings which remain 
after applying the pulse sequences described above are 
those between qubits with different resonant frequencies 
{i 7^ j) and in different chains {m ^ n). These couplings 
only arise when two planes of qubits are selectively re- 
coupled during a two-bit gate. The needed time duration 
for a logic operation may be minimized by maximizing 
the desired qubit-qubit coupling - this is accomplished 
by first letting a be the smallest nuclear distance in the 
crystal. Second, we set = 0, hence putting the gradient 
parallel to the applied field. Third, only nearest-neighbor 
qubit couplings are used (j = J + f ), with longer distance 



couplings accomplished by bit-swapping IQ . After mak- 
ing these choices, we have 5u! = Jwi.i+i — — /xo7^?i/27ra^, 
and the target qubit of frequency Ui in the ?77,th chain 
is coupled to control qubit copies at oj^+i in the other 
chains by the Hamiltonian 
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The decoherence caused by this undesired Hamiltonian is 
minimized if the coupling constants are sufhciently small 
to obey 
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where a (the square root of the second moment 15|) 
is the effective linewidth of the qubit ensemble during 
recoupling, which we demand to be much smaller than 
the splitting 5lo. This inequality is best satisfied by a 
crystal whose nuclei couple only in isolated chains. 

A crystal which approximates this description is fluo- 
rapatite, Ca5F(P04)3, whose ^^F structure is shown in 
Fig. |l|. In fluorapatite, we find u/Slo ~ 1/58, roughly 
six times smaller than in a simple cubic crystal such as 
CaF2. The one-dimensional nuclear structure of fluora- 
patite has been recognized in several NMR experiments 
|l7| . Decoherence timescales in this crystal are analagous 
to the well-known case of CaF2. The T2 time is lim- 
ited by dipolar broadening; if dipolar couplings are per- 
fectly controlled, the timescale for internal decoherence is 
pushed toward the spin-lattice relaxation time constant 
Ti, which is limited by thermal fluctuations of paramag- 
netic impurities (rare-earth substitutions, for example) 
and can easily be several hours for reasonable crystal pu- 
rities and temperatures ||l8| . 

The field gradient is accomplished by a micron-sized 
ferromagnetic parallelopiped. Calculations similar to 
those in Ref. |lj] show that such a magnet, made of dys- 
prosium and placed in a 7 T external field as shown in 
Fig. H, can produce a field gradient of 1.4 T//xm at a 
distance of 2.07 /zm above the magnet, which leads to 
Aw = Awi^i+i = 27r X 19.2 kHz. A 1 /xm by 10 ^m 
fluorapatite crystal above the magnet contains N = 10^ 
equivalent-frequency qubit copies in each xy-plane. The 
field variation along the y-direction is negligible since 
the magnet is much longer than the crystal. The field 
variation along the a;-axis is sufficiently small that all of 
the equivalent-frequency nuclei lie within a bandwidth 
of Aw. We also note that inhomogeneous broadening is 
constantly refocused in this scheme by the narrow-band 
IT pulse sequences during both decoupling and selective 
recoupling. 

The presence of a large magnetic field gradient provides 
a natural means for performing MRFM on a magnetiza- 
tion M^, since this technique is sensitive to the gradient 



force given by F = M^Vi?^ [Q. The crystal is mounted 
on a microfabricated cantilever which oscillates in the 
z-direction, as in Fig. |l|. The long axis of the crystal- 
cantilever heterostructure and magnet are aligned. The 
experiment is performed in high vacuum (< 10~^ torr) 
and at low temperatures. A coil is used to generate the 
RF pulses for logic operations and decoupling sequences; 
it also generates the continuous-wave radiation for read- 
out. An optical fiber-based displacement sensor is used to 
monitor deflection of the cantilever using interferometry. 
Sub-Angstrom oscillations can be detected; larger oscil- 
lations can be damped with active feedback which avoids 
additional broadening while maintaining high sensitivity 

Readout is performed using cyclic adiabatic inversion 
p2[ , which modulates the magnetization of a plane of 
nuclei at a frequency near or on resonance with the can- 
tilever. The spins of resonant frequency Ui are irradiated 
with the RF field B^ — 2Bi cos{u)it— {Q/uJm) cos{uJmt)}, 
where uj„i is the modulation frequency chosen to be near 
the resonance of the cantilever, and f2 is the frequency 
excursion, which should be much smaller than Auj [GSJ. 
Simultaneous detection of signal from multiple planes is 
possible if the different planes to be measured are driven 
at distinct modulation frequencies w^. 

The force resolution for MRFM is limited by thermal 
fluctuations of the cantilever |24|. Force resolutions of 
5.6 X 10^^* N/VHz have been reported for single crystal 
cantilevers at 4 K ||2^. To estimate the force involved 
in a measurement following a quantum computation, we 
consider the case in which initialization is imperfect, so 
that the temperature T is non-zero and an effective pure 
state must be used |g|,0. The usable magnetization is 
estimated H as 



M^ = jh 
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In the high temperature limit, this magnetization ap- 
proximates to {-f^ti^ Bq /2kBT)Nn2~"' showing the expo- 
nential downscaling which plagues solution NMR. Even 
in this limit, with the relatively small N of 10^ and using 
the field gradient described above in a 7 T field at 4 K, 
the force which needs to be measured works out to about 
10~^^n2~" newtons, allowing n ~ 10 qubits. The density 
matrix of such a system would be outside the neighbor- 
hood in which it has been shown that all density matrices 
are separable [g|. As seen in Fig. g, the situation is im- 
proved in larger fields and at lower temperatures, where 
Af ^ improves exponentially in Bq/T. In the realistic lab- 
oratory extreme Bq = 20 T and T = 10 mK, this design 
may scale to n ~ 300 qubits; in this regime entanglement 
is shown to be demonstrable B| . 
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FIG. 3. A plot showing the scalabihty of the present 
scheme. The soUd curve, corresponding to the left axis, 
shows the magnetic field (in T) divided by temperature (in 
K) needed in order for a number of qubits n to be measur- 
able. The dashed lines, corresponding to the right axis, plot 
the number of logic gates (decoherence time T2 divided by the 
pulse sequence cycle time ic) times the length of a "block" of 
the decoupling sequence L against n for several values of T^ . 

A primary limitation to the long-term scalability of 
this scheme is the number of logic gates which can fit 
into the decoherence time when there is a large number 
of qubits. The decoupling/recoupling pulse scheme based 
on Hadamard matrices has a cycle time tc{n) of Lri^/Auj, 
where the parameter L depends on the length of the 
homonuclear decoupling subsequence and only weakly 
on the number of qubits n JMJ. Since single spin ro- 
tations must occur between cycles, the clock speed is set 
by this cycle time. The number of logic gates we may 
perform is approximately tc{n) divided into the remain- 
ing decoherence time Tj". This timescale T2 is due to 
residual homonuclear couplings, cantilever thermal drift, 
and thermal relaxation resulting from magnetic impuri- 
ties; these contributions are limited by the experimen- 
tal performance of the homonuclear decoupling pulse 
sequence, cantilever feedback stabilization, and crystal 
growth technique, respectively. Figure shows the num- 
ber of logic gates possible for several values of T^. The 
bottom dashed trace (Tj = 100 ms) corresponds to the 
case in which homonuclear couplings are perfectly con- 
trolled but cantilever drift is completely unsuppressed, 
while the middle {T§ = 10 s) and top {T§ = 1000 s) 
traces assume suppression by 20 and 40 dB with a nega- 
tive feedback circuit. The limited number of gates at high 
numbers of qubits is a remaining problem which may be 
addressed through error correcting codes Q, more so- 
phisticated pulse sequences, or the use of a more one- 
dimensional system [E6|. 

A remaining challenge in this scheme is the alignment 
of the field gradient, cantilever, and crystal. Such align- 
ment is critical, and although difficult, it is not impossible 
with modern microfabrication technology and feedback 
control using micro-electro-mechanical actuators p7| . 

This proposal features two important advantages over 
schemes using solid-state NMR with impurity dopants 
[P8| . The fabrication difficulty of artificially implanting 
controlled arrays of spins is avoided in our proposal by 



using nuclei that are naturally organized into a crystal 
structure. Also, the use of an ensemble of 10^ spins avoids 
the need for the daunting task of reading and initializing 
single nuclear spins. Indeed, ensemble measurement has 
given solution NMR a developmental head-start against 
other existing proposals for quantum computation, and it 
is hoped that the scheme presented here will carry NMR 
quantum computation onward to the many-qubit regime. 
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